Infrared self-consistent solutions of bispinor QED3 
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Quantum electrodynamics in three dimensions in the bispinor formulation is considered. It 
is shown that the Dyson-Schwinger equations for fermion and boson propagators may be self- 
consistently solved in the infrared domain if on uses the Salam's vertex function. The parameters 
defining the behavior of the propagators are found numerically for different values of coupling con- 
stant and gauge parameter. For weak coupling the approximated analytical solutions are obtained. 
The renormalized gauge boson propagator (transverse part) is shown in the infrared domain to be 
practically gauge independent. 
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I. INTRODUCTION 

QED in three space-time dimensions has become a 
testing laboratory for certain nonperturbative aspects 
of QFT like chiral symmetry breaking [ll-[lol|. bound 
states [IMl or confinement % 14- 18] and for various 
approximation schemes. Due to the dimensionality of the 
coupling constant (as -y/mass) it is a superrenormalizable 
theory so it is free of infinite renormalization ambiguities 
typical for four-dimensional theory. 

It can be formulated in two inequivalent versions: with 
two- and four-component fermions [19, 20]. The proper- 
ties of the theory are different in these two cases. To 
the investigation of both versions much attention has 
been payed over the last twenty years. The work has 
been especially concentrated on nonperturbative solu- 
tions of Dyson-Schwinger (DS) equations with different 
approximations incorporated in the theory as quenched 
approximation, rainbow approximation, l/TV expansion 
and various models of vertex function [ll-I3, @, [13 ■ Par- 
ticularly often the multi-flavor theory with the limit 
TV — ^ oo has been used, since it avoids infrared prob- 
lems, which usually become troublesome in lower num- 
ber of dimensions. The other important point in this 
analysis has been the unpleasant gauge dependence of 
the nonperturbative results which constitute the com- 
mon problem of the approximated studies based on DS 
equations S i, [13 [11 [II [H . 

As it is well known, DS equations constitute an in- 
finite set of relations involving Green's functions which 
form the 'inverted ladder' type structure: the n-point 
functions depend on n -f 1-point ones and so on up to 
infinity. This set cannot be solved without the trunca- 
tion of such a hierarchy. But even such truncation, which 
turns the infinite set of equations into only a couple of 
them, leads to the system which is far from being trivial 
and requires further simplifications. Usually this trunca- 
tion is accomplished by performing certain assumption 
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for the vertex function, which is chosen most often in the 
form satisfying Ward-Takahashi identity. Naturally this 
identity fixes only longitudinal part of the vertex leaving 
the transverse part a subject of further discussion and 
improvements [15, 23-27]. 

In our previous paper |28] we applied a method elabo- 
rated earlier in QED4 |29j to QED3 with two-component 
fermions. It consists on the following five steps: 

1. a certain infrared form of two basic propagators, 
S[p) and D^''{k) (suggested by perturbative calcu- 
lations or other methods) as dependent on a couple 
of unknown parameters is assumed, 

2. the fermion propagator is represented in the spec- 
tral form with one known spectral density p{M), 
which is in general possible in the infrared domain, 

3. Salam's form of the vertex function [13], together 
with S{p) and D^'^{k), is substituted into the first 
two of the set of DS equations, 

4. from these two equations the set of self-consistent 
equations for parameters is derived, 

5. the obtained equations are solved numerically or 
analytically. 

This method proved to be relatively effective both in 
QED4, where all parameters were correctly found with- 
out the necessity of infinite renormalization, and in spinor 
version of QED3, where obtained results stay in general 
agreement with other works. In the present paper we 
would like to extend its application to four-component 
QED3, without Chern-Simons term, i.e. when gauge 
bosons ('photons') remain massless in spite of interac- 
tion. 

The common effect of various simplifications of DS 
equations is the gauge dependence of the results (even 
of the physical observables). The sources of this undesir- 
able behavior are the approximations made to the prop- 
agators and to the vertex. It seems therefore valuable 
to test various possible approaches with regard to that 
particular feature and much work has already been done 
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in this direction (see the references above). This point 
Ues in the scope of interest of the present work too. 

This paper is organized as follows. In the next section 
we define the model itself and give the resulting set of DS 
equations. In section Hill the infrared Green's function in 
question are formulated up to several unknown parame- 
ters. In section HVl we substitute these Green's functions 
into DS equations and obtain the set of relations for the 
introduced parameters. In the last section we present 
analytical and numerical results and some conclusions. 



II. FORMULATION OF THE MODEL. 

The model is defined through following Lagrangian 
density : 

C{x) = (17^9^ - mo - eoj^A^ix)) (1) 

where A is the gauge parameter. The quantities mo and 
eo denote here the bare fermion mass and the bare cou- 
pling constant respectively. As mentioned in the Intro- 
duction, the latter for D = 2 -I- 1 is a quantity with the 
dimensionality of y^inass. That means that in the quan- 
tum theory higher terms of perturbation expansions have 
better ultraviolet momentum dependence in loop integra- 
tions and the model is superrenormalizable. 

As already told, in the present paper we deal with four- 
component fermion field, choosing the following represen- 
tation for gamma matrices used also in four-dimensional 
QED: 
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There are two other gamma matrices, which anticom- 
mute with all above, and which can serve for defining 
chiral transformations; 



7^ = 





M 
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(3) 



The fermion mass term, that was chosen in ([1]) in the 
form mo^'^', breaks the chiral symmetry defined by any 
of the matrices (jS]). There is, however, the possibility of 
the other choice for the mass term 



rriQ'^T'^. where 



[7^7^] 



1 
-1 



(4) 



This term does not break chiral symmetry, since matrix 
7°r commutes with both 7"^ and 7^, but it does violate 
parity symmetry (contrary to the choice made in ([1])). 
For matrices ([2]) we have ordinary relations: 

{7^7''} = 5'^% tr7^ = 0, Tr[7^'']=V^ 

Tr [7^7"^7''7l = 4 {g'^-'g'"^ - g'^Pg^" + g^'^g'^") , (5) 



where we choose for the metric tensor: 
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The trace of the product of an odd number of gamma 
matrices equals zero as in four dimensions (this was not 
the case in the spinor representation, where it was pro- 
portional to the antisymmetric tensor e^"'^). Addition- 
ally we have the identity: 
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Using the Lagrangian density ([T]), one can derive in 
the standard way - for instance through Feynman path 
integral - the Dyson- Schwinger equations for propaga- 
tors (sij . For boson propagator we obtain the relation 
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sip)rp{p,p-k)s{p-k)D^^k) 



(8) 



which may be given the graphical form shown in Fig- 
ure [TJ The propagator D^'^{k) on the right hand side is 
not integrated over three-momenta and therefore it can 
fully be represented through fermion functions (and the 
vertex) . 
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FIG. 1: Dyson-Schwinger equation for the gauge boson prop- 
agator D'^"{k). Light lines represent free propagators and 
heavy ones dressed propagators. The full circle stands for the 
full fermion-boson vertex. 

The DS equation for the fermion propagator does not 
allow for such a separation, and has the form 
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S{p + k)T''{p + k,p)S{p)D^,{k) 
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FIG. 2: Dyson-Schwinger equation for the fermion propaga- 
tor. As in Figure [Tl heavy lines stand for full functions and 
light for free ones. 

Its graphical representation is shown in Figure [H 
While introducing the fermion self-energy S(p), it can 
be rewritten in a simpler manner 



{^-mo)S{p)^l + nP)S{p) 



(10) 
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where 



(11) 
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s{p + k)r''{p + k,p)s{p)D^,{k) 



Higher equations are not considered in our present 
approach. The vertex function F'^ is not taken from 
DS equations, but postulated in the form proposed by 
Salam [30| and used afterwards in the so cahed 'gauge 
technique' |24, 25, 32, 33]. The eventual extension of our 
method on higher Green's functions will be considered 
elsewhere. 



III. INFRARED GREEN'S FUNCTIONS 

The free propagator of massless vector boson has the 
standard form: 



-g 



1 kf'k" 



A (fc2 



(12) 



As is well known, the interaction with fermions does 
not change the longitudinal part of Z?'"^(fc). This is guar- 
anteed by the following Ward-Takahashi (WT) identity 

A^i?^''(fc) = fc^i^(")'-(fc)--i^. (13) 

Consequently the nonperturbative propagator may be 
then written as 



D'^'ik) 



d{k^ 
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(14) 



with certain unknown function d{k'^). In the infrared do- 
main we assume this function to have a Taylor expansion 
starting from fc^, since for bispinor fermions the boson 
field remains massless J^] and we expect d{Q) — 0. This 
is not the case for massive two-component fermions, for 
which the topological boson mass is generated. 

Taking the first two terms of this expansion, we have 



(15) 



d{k') = fc^ 1 + — 



where k^, together with the renormalization constant Z3, 
will be determined from consistency conditions. From 
the unitarity we expect the value of to satisfy the re- 
striction < Z3 < 1 [s^l- This expectation was actually 
confirmed in our previous works on spinor QED3 [28j and 

QED4 m. 

With the above assumptions the inverse of D^^'^ may 
be written as 



(16) 



1 + — (-fc^g^'^ + fc^r) - A/c^fc" 



The free fermion propagator has the usual form: 



i - Too 



Due to the masslessness of the gauge boson and the ab- 
sence of the mass gap for emission of soft photons, the 
pole at mass too should, in the full propagator, turn 



into a branch point at p 



where to is a phys- 



ical mass. This is expected also on the basis of gen- 
eral considerations on the analytical structure of fermion 
propagator [l3, ll^ and of confinement, which prohibits 
S{p) from having a simple pole. The numerical results, 
performed in 1 /N expansion with bare vertex and in eu- 
clidean space, suggesting the existence of complex sin- 
gularities rather than real may as well constitute the 
effect of the coarse approximations made while solving 
DS equations: the approximations that are inevitable in 
any nonperturbative approach. Such singularities may 
be a signal of confinement, but are not prerequisite. In 
the Schwinger Model, massless electrodynamics in two 
space-time dimensions, which exhibits confinement of 
fermions, the real singularity in the infrared domain in 
the fermion pro pag ator has been found (it has the form 

of i/{~p'f/^)mM- 

Our infrared assumption for S{p) is then: 



S{p)^ 
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(18) 



We will see later that it will turn out to be self-consistent. 

The values of the exponent /3 and of mass renormal- 
ization constant 8m = to — toq will be established by the 
requirement of consistency. The presence of fi in denom- 
inator improves the ultraviolet behavior of loop integrals 
(we assume that < /3 < 1, to be verified a posteriori). 

For our further purposes S{p) has to be written in the 
spectral form 



S{p) = J dMp{M) 

with the one spectral density 
sin(7r/3) 
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(19) 



p{M) 



TT (Af -m)(M2/m2 - 1)^ 
X [e(M - TO,) - e(-M - to)] , 



(20) 



sufficient to define the infrared behavior of the propaga- 
tor. 6 is here the Heaviside step function. The particular 
form of p{M) was given in our previous works [H, H^. 

Using p{M), one can write the vertex function, which 
we need to put into DS equations ^ and ©. As men- 
tioned in the Introduction we use the Salam's vertex [s^] 
(with slight and obvious modification resulting from the 
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form (|19[) of S{p), which in our case contains two terms): 

sip + k) r^ip + k,p)sip) 

dMp{M) 



— m p — m 
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(21) 



It automatically guarantees the compliance with the WT 
identity: 

k^S{p + fc)r^(p + k,p)S{p) = S{p) - S{p + k) . (22) 



IV. SELF-CONSISTENT INFRARED 
EQUATIONS 

A. Gauge boson propagator 

Inserting and (PT|) into the right hand side of 
we obtain 
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M + ie 



m + ie 
(23) 



where we have rewritten this equation for inverse propa- 
gator, which is easier to handle, because in this case D^"^ 
decouples from other functions. 

The vacuum polarization tensor may be defined as the 
spectral integral 



n^^(fc) = J dMp{M) {W^{k) - nZ{k)) . (24) 
where n^J^(fc) denotes the usual perturbative tensor: 



n^^'^(fc) = 



(25) 



(Pp 



(27r)3 ^ — m + ie ^ — ^ — m + is 



It may be evaluated in the standard way with the use 
of identities ([5]), and for instance by performing Wicks 
rotation and introducing Feynman parameters. Passing 
back to Minkowski space, we have 



(26) 



fl(_fc2gA'!^ + k^^k") [ da;^ . 

TT J (tit? — k'^x{l ~ x)y/'^ 



For k^ < 4m^ the x integral is well defined. The 
trans versality of n^J'(fc) is a favorable consequence of us- 
ing the vertex function in the form (|2ip. satisfying the 
WT identity (|22l) . Performing the spectral integral over 



M in f24)) with the use of (|20|) . similarly as it was done 
in 1281, we obtain: 



m^ik) 



eW + 1/2) 
7r3/2r(/3 l)m 



(27) 
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The integral over Feynman parameter x leads to the 
hypergeometric (Gauss) function 



dx ■ 



x{l ~ x) 



(l-2/x(l-x))^+i/2 

i2i^i(2,/3+l/2;5/2;y/4), (28) 
6 

and the DS equation ([25)) may be given the form 
D'^ikY" = -Xk^k" + i^k^g^"" + k'^k'') X (29) 
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67r3/2r(/3-f l)m 



2Fi(2,/3 + l/2;5/2;fcV4TO2) 



After the substitution of the expression (IT5|) for the left 
hand side and cancellation of the tensor structures, we 
are left with the scalar equation, for which we require the 
adjustment the first two terms of the Taylor expansion 
in fc2. In that way we get two equations for and k^: 



r,-! _ . , elnp + 1/2) 
3 + 67r3/2TOr(/3 -f 1) 



egr(/3 + 3/2) 
k2 307r3/2m3r(/3 -I- 1) 



(30) 
(31) 



where we expanded the Gauss function for small momen- 
tum according to the formula: 

2Fi{a,h;c]z) « 1 + —z + 0{z^) , 
c 

and r is the Euler function. These are two of the set of 
four equations for parameters of the model, to be solved 
in section Ivl 



B. Fermion propagator 

To get other two equations for parameters, we put (|19l) . 
(|2T|) and ^ together with ([15]) into the DS equation 
We obtain, after some simplifications 



(|(- mo)S'(p) = 
l + [np)S{p)]A + [np)S{p)]B 



(32) 



[np)S{p)]c 



where we divided fermion self-energy into pieces com- 
ing from different tensor structures in D^'^^^^: g^^ and 
k^'k'^ from the transverse part and again from the gauge- 
dependent longitudinal part. After small rearrangement 
they are: 
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[ms{p)]A = 

iZ3K^el J dMp{M) 

[np)s{p)]B = 

iZ:iK^el j dMp{M) 
[nP)S{p)]c = 



(27r)3 



1 



1 



^ + ^ — m + ie (fc2 + ie){k'^ — + ie){-^ — m + ie) 
+ ^ -m + ie){k'^ + ie)'^{k^ - + ie) ~^ 



d^k 
^7 



— - (m ^ M) 



(33) 
(34) 
(35) 



1 



+ ji — m + is){k'^ + ie)"^ 



(m Af) 



A comment should be made here. To avoid technical 
difficulty while performing Wick's rotation in the above 
momentum integrals, we analytically continued the value 
of K to imaginary values on the upper half plane of com- 
plex At. Now the deformation of the integration contour 
as is required by passing into euclidean space, is not dis- 
turbed by the inappropriate location of poles since all sin- 
gularities have the 'Feynman' position. After performing 



the integrals, we will come back to real values of k. The 
procedure was discussed in [2^. With this trick each 
of the above momentum integrals can be performed in 
an ordinary way known from perturbation theory. It can 
easily be seen that all integrals are finite without the need 
of any regularization. Omitting details of this standard 
calculation, we find (in Minkowski space) 



I A = iCQ 



d^k 
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(27r)3 ' ^ + /(f - m + ie^'^ (fc^ -|- ie){k'^ - + ie) 
1 



(36) 



Ib = ieo 



dx(3m— ^(1 — x)) 
d^k „ 



1 



1 



{m'^x — p^x(l — x))^^^ [m?x — p^x{l ~ x) + k^{1 — x))^^-^ 

1 



(27r)3 ' {^+ ^-m + ie){k'^ + ie)'^{k'^ - + ie) 8ttk^ 
1 1 



dx 



1-x 



+ m) 



(37) 



{m'^x — p^x{l — x))^^^ [m^x — p^x{l — x) + k^(1 — x))^/^ 

1 



+ to) x{1 ~ x) 



2 (m^x — p^x(l — .t))3/2 



Ic = 



d^k 



1 



= dx 



(27r)3 ' (|< + - m + i£)(fc2 + ie)2 Stt 
+ to) x(1 — x) 



{rii?x ~ p^xll — x))^/'^ 



2 (to^x — p2x(l — a;))3/2 



(38) 



r 



Now the contributions to Yi{p)S{p) may be written as All the above spectral integrals can been found similarly 

as in |2a, therefore we omit the technicalities. 



[np)s{p)]A = 

ZgK^y" dMp{M) 



Ia- 
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m + ie 



(to -> M) 



(39) 



(40) 



[S(p)5(p)]c = J I dMp{M) [Ic - (to ^ A/)] . (41) 



[np)S{p)]B = 
ZgK^ / dMp{M) [Ib - (to ^ M)] , 



In the infrared domain, when — > , the left hand 
side of the DS equation (l32t , after substituting ([T8| , con- 
tains two types of singular terms (the only terms that are 
important): 



i'^^Sm{^ + m) TO 



2/3 



(to^ — p^'J^i + l ^jji2 _ p2'jl3 



(42) 
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Therefore, to avoid lengthy expressions we wih not give 
the results for (gH]), and dH]) in their full complexity, 
because it is sufficient for our goal to only pick out from 
them the identical singular terms. Wc find: 



— Ii(m^, K^) 



+ m) 



(f -p7m2)/3+i 



■^{^ + m) 
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+ m) 1- P 
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1 \ 1 



(43) 
(44) 



[np)S{p)]c 



2m3 13 4:(3mJ 

Cq [^(^ + m) f 

+ m) -1 



m3 (l-p7m2)/3+i 
1 \ 1 



2m^ 



(45) 



P A/Sm J (1 -p2/m2)'3_ 

where w refers to diverging terms, when — m?. The 
functions Ii and I2 have the following form: 

^-±1 , (46) 

{rn^x^ + ^2(1 — a;))-'-/2 



Ii(m^,K^) = j dx 



SOLUTIONS AND CONCLUSIONS 



The four equations we have obtained, i.e (pO)) . (|3T1) . 
(HHl) and (15(11) are sufficient to determine all parameters. 
We rewrite them with the use of renormalized quantities: 

fermion mass m = mo + Sm, gauge coupling constant 

1/2 

e = Z3 eo and gauge parameter Xr — Z3X. Besides, it 

2 

is useful to introduce a dimensionless parameter C — ir — . 

^ ^ 47rm 

After executing the parametric integrals in Xi and Z3 and 
performing the reverse analytical continuation in k, we 
get 



6m 
m 



4to2 



^-1 



4/3 V Ak 



4to2 



ln(4mVK^ + 1) 



4m^ 



4m2 + k2 



ln(4mVK^ + 1) 



2m2 

^ 2C r(/3 + l/2) 
3 y^T{l3 + 1) 
15 v/H-(/3 + l) 
2C r(/3 + 3/2) 



,(52) 

(53) 
(54) 
(55) 



This set of equations for Sm/m, /3, Z3 and jw? may 



T2(m^, K^) = dx 



x'^ix + 2) 



(m2a;2 + k2(1 — a;))'^/2 



(47) 



If the solution is to be self-consistent close to the 
fermion mass shell, the divergent terms on both sides 
of DS equation (ITU|) must be identical. Equating them, 
and making use of the fact that up to finite terms we 
have 

^(|( + m) 2to2 



{1-pym^y [l-p^/m^Y ' 
i + m) rn(^ + m) 



(48) 



m 



(l-p2/^2)/^+l (l_p2/^2)/3 + l (i_p2/^2)/3' 

we derive the following two relations for unknown param- 
eters 5m and /3: 



Stt 
Stt 



-mli(m2, ^2) -f 

AZ3 



4pm \AZ3 



(49) 



(50) 



The function 13(771 , k ) is defined as follows 
— (Xi (777^, K^) — 777^X2(777^, K^)] 



(1 -a:)(x + 2) 
(7772a;2 + k2(1 — a;))'^/2 



(51) 




FIG. 3: The dependence of k in units of m (upper plot) and 
of power P (lower plot) on the parameter C,. The dashed 
line corresponds to approximated solutions defined by equa- 
tions (f55)) -((5 ^ . The gauge parameter is chosen as Ah = 0.3. 
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be solved numerically for certain chosen values of renor- 
malized gauge parameter Xji and the results plot as func- 
tions of parameter (. Let us consider the case of weak 
coupling, when is small. Since the applicability of our 
method requires < /? < 1, then from equation ((55)) 
we deduce that /m? should be large. It may be jus- 
tified by the elementary estimation given below. First 
we rewrite following expression in terms of beta function 
B{x, y) and its integral representation: 



^T{f3 + 1) _ 1 
2r(/3-h3/2) 2 



B 



(/3 + l,l/2) = i / t-^l^(l~tfdt 



and next use the inequalities valid for < /3 < 1: 

1 1 



1 1 



6m 





Now, from ((55)) it becomes obvious that for small ^ the 
left hand side must be large. This is in agreement with 
our expectations concerning the subsequent terms in the 
Taylor expansion ((T5)) . But for ^ the first term in 
square brackets on the right hand side of ((53p dominates 
over all other, which may be easily verified. The expected 
positivity of (3 requires then considering only gauges for 
which < Xr < 1/2. This is not surprising, since /? is 
a strongly gauge dependent quantity. For our numerical 
calculations we have then chosen the values of the gauge 
parameters from that range. 

In figure [3] we show the dependence of parameters k 
and P on C for exemplary value of Xr = 0.3 (solid lines). 
In figure |4] we similarly plot the dependence of two other 
parameters: fermion mass renormalization constant Sm 
and gauge field renormalization constant Z3. It is nice 
to observe that we have < 6m < m and < < 1, as 
we expected. The dashed lines in these figures represent 
approximate solutions of the set ((5^ - ((55)) . as described 
below. 

In the weak coupling regime, assuming that we do not 
choose the gauge parameter approaching zero (but 
still to be less than 1/2), the equations may be given the 
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FIG. 4: The dependence of mass renormalization Sm in units 
of m (upper plot) and charge renormalization constant Zs 
(lower plot) on the parameter ^. The dashed line corresponds 
to approximated solutions defined by equations (|56p - (|59|) . On 
the upper plot the dashed line is almost identical with the 
solid one and therefore it is not visible. The gauge parameter 
is chosen as = 0.3. 



FIG. 5: The comparison of the behavior of k (upper plot) 
and parameter f5 (lower plot) for different values of gauge 
parameter: dotted line - A_r — 0.1, dashed line - Xr = 0.2, 
solid line - A_r = 0.3, mixed line - A^ = 0.5 (Yennie gauge). 
The plots of K for various gauges follow almost the same curve. 
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5m 



0.1 0.2 0.3 



The terms / in the first two equations arise from 
the expansion of Ti and I3 for large , but in fact they 
may be omitted, since they are of order as results 
from the last equation. As already told, the solutions of 
these simplified equations are drawn in figures [3] and |4] as 
dashed lines. On the plots for k and 5m these curves are 
not visible since they are almost identical with the full 
solutions, but also on the other two they do not deviate 
from the 'exact' results too much. 

From (1571) we see that the exponent /3 in this approxi- 
mation may be written as 




(60) 



0.2 0.3 



FIG. 6: The comparison of the behavior of fermion mass 
renormalization 5m (upper plot) and gauge field renormal- 
ization constant (lower plot) for different vales of gauge 
parameter: dotted line - \r = 0.1, dashed line - \r — 0.2, 
solid line - \r = 0.3, mixed line - Xr = 0.5 (Yennie gauge). 



following approximated form; 
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(56) 
(57) 
(58) 
(59) 



which means that Xr = 1/2 is a kind of Yennie 
gauge. 

It is interesting to observe, if and how the values of 
the parameters depend on the gauge. This is shown in 
figures [S] and [S] which are performed for the following 
values from the range [0,1/2]: = 0.1, 0.2, 0.3, 0.5. 
The Landau gauge cannot be used because it would lead 
to negative value of f3 . 

The particular stress deserves the observation that the 
gauge dependence of the parameter k is extremely weak. 
This means that the renormalized gauge boson propaga- 
tor (strictly speaking its transverse part) is practically 
gauge independent, as it should be. Please note that 
in other approximation schemes applied to DS equations 
in QED3 one obtains the gauge dependent value of the 
polarization scalar [l6j . 

The dependence of on is relatively weak for small 
C too. These results are worth noticing, since the full 
gauge independence should appear in the exact theory, 
and one ought not to expect too much from the approxi- 
mated model, where infrared forms of Green's functions 
are postulated in a simple form. Gauge dependence of the 
'physical' fermion mass obtained in our work is relatively 
strong but it is a common feature of nonperturbative cal- 
culations in this model [1, 0, [2l| . 
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